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Abstract 

We study the properties of rotating neutron stars within a generahzed chiral SU(3)-flavor model. 



discussed. We calculate the Kepler frequency and moments of inertia of the neutron star sequences. 



■ An estimate for the braking index of the associated pulsars is given. 
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I. INTRODUCTION 



Neutron stars belong to the most important physical systems where strong interactions 
are "tested" at extreme dense conditions. In contrast to ultrarelativistic heavy-ion collisions 
where also high baryonic densities might be obtained by choosing the appropriate bom- 
barding energy, neutron star matter exists at practically zero temperature - at least on a 
hadronic scale. Thus the understanding of neutron stars is the complementary task of what 
lattice gauge calculations can achieve for high-temperature, zero-density matter. In addition 
due to the fact that charge neutrality drives the stellar matter away from isospin-symmetric 
nuclear matter the study of neutron stars can give important clues for understanding the 
isospin-dependence of nuclear forces. A related regime of strong interactions can be exper- 
imentally investigated using secondary beams of nuclei going towards the neutron drip-line 
of nuclei. 

Compared to the harder task of determining masses or radii the most easily observable 
property of a neutron star is its rotational frequency through its pulsar radio signals emitted 
from the magnetic poles of the star[0]. Additional information can be obtained by the slow 
change of the pulsar frequency of an isolated neutron star with time signalling the slowing 
down of the rotation due to radiative loss of angular momentum. So far observational data 
in this respect are still scarce ||, ^ 0. 

Various hadronic models have been used to describe the structure of neutron stars. Going 
to central densities of a few times nuclear matter saturation density the inclusion of hyperon 
degrees become essential [0]. As it always will be difficult to get clear and unambiguous in- 
formation on the stellar interior from the relatively limited observational data it is important 
to relate star properties with other observables like neutron distribution radii in neutron-rich 
nuclei or multi-fragmentation in heavy-ion collisions with different N/Z nuclei. It is there- 
fore very useful to implement a model with a broader applicability including nuclear matter, 
finite nuclei and chiral symmetry restoration to be able to simultaneously study different 
observables within a unified approach using the same set of parameters for all calculations. 
In this paper we follow this direction by adopting a hadronic model incorporating the un- 
derlying SU L{'i)xSU r{?>) structure of QCD, which we used in earlier studies of non-rotating 
neutron stars, nuclear matter and finite nuclei/hypernuclei calculations [0, ^J. 

Here we extend the static calculation to study the properties of rotating neutron (hyper-) 
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stars within the same approach. The outhne of the article is as follows. First we give a brief 
description of the ingredients of the hadronic model that we implement in our calculations 
(a much more detailed discussion can be found in |]10|)- After introducing the Tolman- 
Oppenheimer-Volkov (TOV) equations for the static star we present the extension to rotating 
stars following the treatment of Hartle and Thorne In the third section we present some 
general features of the equation of state and its strange matter contributions and then look 
at the composition and gross properties of the rotating stars. We conclude with an outlook 
of possible extensions of the current approach. 

II. THE CHIRAL MODEL 

The chiral hadronic SU(3) lagrangian has the following basic structure 

= ^kin + -^BM + -^BV + -^vec + -^^0 + -^SB + -^lep ) (l) 

consisting of interaction terms between baryons and spin-0 (BM) and spin-1 (BV) mesons 

i 

Aec = Iml^co' + lml4^' + \t,my + gl{u' + 20^ + Qu'p' + p') (2) 



summing over the baryonic octet, (p,n,S ^^), as well as interactions between the 

scalar mesons 



-hx --,X ln^ + -ln^— . (3) 



An explicit symmetry breaking term mimics the QCD effect of non-zero current quark masses 
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Finally, there are leptonic contributions from electrons and muons 



(4) 



Aep = 'ipiii^i^d'' - mi]'^i . (5) 

l=e,fj, 

£km in (1) contains the kinetic energy terms of the hadrons. The general model incorpo- 
rates the complete lowest baryon and meson multiplets. In the previous formulae we only 



considered the degrees of freedom relevant for a neutron star calculation, the scalar field 
a and its ss counterpart C (which can be identified with the observed /o particle), as well 
as the uj,p and vector mesons. £vec generates the masses of the spin-1 mesons through 
the interactions with spin-0 mesons. The scalar interactions Cq induce spontaneous chiral 
symmetry breaking. Another scalar, isoscalar field incorporated in the model, the dilaton 
X simulates the breaking of the QCD scale symmetry and can be identified with the gluon 
condensate|lT3|, for a longer discussion see [|lOl)- The effective masses m*(o", Q = Qa o- + C 



of the baryons are generated through their coupling to the scalar fields, which attain non- 
zero vacuum expectation values due to the self- interactions C, [|10]. As the field strengths 
change in nuclear and neutron matter, the effective masses of the baryon octet are shifted, 
too. The grand canonical thermodynamic potential of the system at zero temperature can 
be written as 

n/v = - £„ - Css - V,. -E^.l'f' l^m - ..-I - 5 S ^ / yi^- (6) 

The chemical potentials for the baryons read 



with bi, Qi, and kp,! being the baryon number, charge, and fermi momentum of the ith 
species. The energy density and pressure are given by 

e = n/V + J2 fikPk (8) 

k=i,l 

p = -Vt/V. (9) 

By extremizing VL/V one obtains self-consistent equations for the meson fields in conjunction 
with demanding vanishing total charge and (3 equilibrium as the time scales involved allow 
for weak interactions to take place. 



III. STELLAR EQUATIONS 



In the case of the static neutron star we integrate the TOV equations|14 



dP{r) 

dr 
dm{r) 

dr 



{e+P){4Trr^P+m) 
r2{l-y{r)) 

: 47rr^e(r) , 



(10) 
(11) 
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where P and e are the pressure and energy density of the nuclear matter. m(r) is the 
gravitational mass of the star inside of the radius r and Y{r) = 2m{r)/r. Starting with 
some initial values for P and e, P{r = 0) = Pc and e(r = 0) = ec(Pc); with m(0) = one 
can integrate equations (|TD|) and ( [TI| ) up to the point of vanishing pressure, which defines 



the radius R of the star. Here the specific model of nuclear matter enters via the equation 
of state e(p). For the low-density tail of the EOS we use a standard parametrization for the 
crust of the star as discussed in p. 

Due to the rotation the star in general is deformed. Assuming axial symmetry the metric 
of the rotating star can be written as: 

ds^ = exp{2u)df - exp{2\)dr'^ - exp(2/x) (dO"^ + sin^ 9{d(j) - udtf) . (12) 

The metric functions /i, A depend on the radius r and the polar angle 9. If we set 
C(j = we recover the standard Schwarzschild metric for the static star with the choice 
exp(2/i) = , exp(2A) = 1/(1 — y(r)). Because of the dragging of the local inertial 
frame uj{r) depends on the radius (but not to lowest order on the polar angle, see [Q) and is 
not equal to the rotational frequency VL of the star. Following |]12| and defining the difference 
of the frequency of the local frame and Vt a.s uo = — uo demanding equilibrium between 
pressure, centrifugal and rotational forces in the local frame the differential equation for the 
metric functions follows [0 



d_ 

dr 



., .duj 4 ., , , , 
j(r) — + -j[r)uj[r) 
dr r 



(13) 



with j = exp(— z/ — A). The rotation of the star generates a modification of the metric as 
shown in ([T2|). Using a multipole expansion of the metric functions, one obtains corrections 
that up to quadrupole order read 

exp(2z/) = exp(20) (1 + 2 {ho + /isPa (cos^))) (14) 

exp(2/x) = r^l + 2{v2- h2P2 (cos 6))) (15) 
fow ^ , 2mo + m2P2 (cos^)\ 

where the static metric component is given by 



The monopole coefficients ho and mo obey the differential equations {r < R): 



drriQ 
dr 

dpo 
dr 



de 



47rr^T7T7(e + P)po + ~3'^ 
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121 -Y \dr 
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^dr\l-Y 



(18) 
(19) 
(20) 



mo{R) 



at the star's surface as can be 



The value is fixed by requiring that hoi^R) = r(i-y{r)) 
inferred from the large-r asymptotics for /io(r) ||15[. Here mo and po are the monopole mass 
and pressure perturbation due to the rotation of the star. The first terms where the actual 
deformation of the star's shape show up are the quadrupole corrections h2, m2 and V2- They 



can be determined by solving |T^, [15 
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(22) 



Integrating the equations and ensuring the correct asymptotics of vanishing corrections 
at r = and r = cxd by adding the appropriate solution of the homogeneous equations 
one obtains the quadrupole corrections to the metric. The mass and pressure quadrupole 
corrections follow as [|I^] 

2 



m2 



P2 



r(l-r) 
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'4) 



,df 



-[ruj) 
3 dr 



ho 



-h2 - — exp (-20) uj"^ 



(23) 
(24) 



After determining the metrical components we can then calculate radii, masses and moments 
of inertia of the rotating solutions. 



IV. NUMERICAL RESULTS 



We use the same parameter set (called "CI" in [0, |T^) that was also used in calculations 
of excited nuclear matter and finite nuclei. Figure 1 shows the result of a calculation of 
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FIG. 1: Strangeness content of nuclear matter as function of charge per baryon and density. The 
full line shows the onset of strangeness in the system, the following lines mark the strangeness 
fraction fs = 0.1,0.2, and 0.3 respectively, where fs is defined as number of strange quarks per 
baryon 



nuclear matter varying the average charge of the system between symmetric nuclear matter 
{Q/A = 0.5) to neutral hadronic matter {Q = 0). One can see how the onset of hypermatter 
is shifted from neutral matter at a density of ~ 2po to twice that value for symmetric nuclear 
matter. Note that here no leptons are taken into account that shift the total baryonic charge 
to positive values in the case of a neutron star. There, A and S states are populated at 
densities higher than 2.6 po 0- Looking again only at the baryonic contributions Figure 2 
shows the the energy per baryon E/ A — rriN SiS a. function of density for total isospin per 
baryon (3rd component) zero and -0.5. The isospin curve exhibits the nuclear matter 
saturation properties. The thick solid line shows the result of the baryonic energy inside 
the neutron star (including leptons) as calculated in the chiral model. One can see that for 
higher densities the equation of state resembles more the isospin nuclear matter results. 
In fact, in the extreme high-density limit assuming similar population of the baryonic octet 
the total isospin tends towards zero again. 

The maximum frequency Qk and the minimum rotational periods Pg = 2-71 /Q^ of a 



rotating neutron star have been determined by integrating the self-consistent equation [15 



u = v{n) exp [u{n) - fi{n)] 



(25) 
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FIG. 2: Binding energy of nuclear matter as function of baryon density. The figure shows isospin 
symmetric matter (dashed line) and neutral baryonic matter including strange baryons (solid line) . 




FIG. 3: Minimum rotational periods Pk for stars with different central energy density. The upper 

two curves are results for the chiral model (full line) and the TMl relativistic mean field model. 
The lower curves are the corresponding Newtonian estimates. 

solving for Q at the equator of the star, where the mass shedding occurs first, for different 
central energy densities and total masses of the star. The resulting values are shown in 
Figure 3. Compared to the Newtonian value of the Kepler period (dashed lines) there is 
an increase in of about 25% due to the general relativity corrections. This result is 
in accordance with other calculations of different neutron star models, which show similar 
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FIG. 4: Mass of different neutron star as function of rotational frequency. The sharply rising 
dotted curve show a result for the TMl parameterization. 



results |T^. In the plot also results for a relativistic mean-field parametrization (TMl [pTl ) 
are shown, which was extended to include hyperons (see @). The difference of the Kepler 
frequency for both parameterizations is small. Note that the maximum gravitational mass 
of the neutron star that was about 1.64 Mq in case of the static solution is now increased 
to 1.94 solar masses due to the additional rotational contribution. Looking at the frequency 
dependence by spinning up the neutron star one gets the increasing neutron star mass 
as shown in Figure 4. One TMl result for a star with M(i7=0) = 1.48Mq is shown for 
comparison. In this case one can see a distinctly sharper rise of the mass of the star with 
rotational frequency. For the chiral case one can observe the increase of the mass by up to 
30% at the mass-shedding frequency. 

Similarly looking at the monopole corrections one can calculate the increase of the radius 
of different star solutions due to rotation averaged over the star's surface. Figure 5 shows 
that there is an increase of the radius by up to 10% close to Qk- Again, the TMl model 
shows a more drastic change with rotation. From calculating the quadrupole coefficients 
Eqs. (|16D one can determine the excentricity eex of the star defined via the ratio of the 
surface radius at the equator over the pole radius: 



-^pole/ -^equator (2^) 
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FIG. 5: Radius of different neutron star (cf. Fig. 4) as function of rotational frequency. 




FIG. 6: Excentricity of the deformed rotating neutron star as function of rotational frequency. 

Figure 6 shows that there is a substantial deformation of the star. This has also an 
impact on the internal structure of the star depending on the polar angle. This is shown 
in Figure 7. Here the onset of the various hyperons contained in the star S^,A and S— 
are shown for a typical static star with a Mass of 1.475 (vertical lines) compared to the 
same solution corrected for rotational effects at its Kepler frequency. 
Integrating the inertia of the star by using the general expression 



J = 27r y" dOdrexp {X + fx + u + ip} {e + P) (exp [2u - 2^] - cu^) ^ (27) 
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FIG. 7: Inner structure of the neutron star. The onset of hyperon population in the interior of the 
neutron star as function of depth and polar angle is shown for A (long-dashed lines), (dashed 
lines) and particles (dotted lines) is shown. The full line represents the star's surface. For 
comparison the thresholds and the radius for the corresponding static star are shown (straight 
lines) . 

and again taking into account the corrections to the metric up to order one gets the 
moment of inertia of the star as shown in Fig. 8. In this case the mass of the non- rotating 
star is 1.48 solar masses. In this figure we have kept the baryonic number of the star 
constant at Mb = 1.22 ■ 10^'' ~ {MB)soiar and changed the rotational frequency. The curve 
therefore shows the change in the moment of inertia during the spin-down of a pulsar, the 
energy-loss of the star generated by electromagnetic and gravitational radiation. Assuming 
a simple power-law of the energy loss AE ~ —E^ where Uo = 3 and rio = 5 in the case of 
magnetic and gravitational radiation, respectively. With this assumption one can determine 
the frequency-dependent braking index n{Q) = rio — An{Q) with 



An 



(28) 



21 + rn 

where the prime denotes the derivative with respect to frequency. Figure 9 shows the shift 
of the braking index. There is a substantial reduction in the case of fast rotating stars that 
has to be taken into account for an analysis of the time evolution of the star. However, one 
cannot see significant structures as for instance a back-bending shape as seen in moments of 
inertia of atomic nuclei and also suggested for neutron stars |0. This is not very surprising, 
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FIG. 8: Moment of inertia of the rotating star as function of frequency. A fixed value of Mb is 
assumed. The result of the chiral model (solid line), the TMl calculation (dashed lines) and in 
comparison the inertia of the star neglecting the deformation (dotted line) in the case of the chiral 
model are shown. 




FIG. 9: Deviation of the braking index from an assumed power law. In the case of magnetic 
radiation the index is reduced by a factor of up to 1/3. 

since the chiral model does not show a phase transition for low temperatures and at high 
densities (see, however, the discussion in the concluding section). Note that in this case at 
maximum rotation the star starts with a central density of 2.3 times nuclear matter density, 
which is below the threshold of the occurrence of hyperons. Below a frequency fl ~ 5600/s 
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(P ~ 1.1ms) hypermatter becomes energetically more favorable. 



V. CONCLUSION 

We have used a chiral SU(3) model for the calculation of rotating neutron stars. The 
model parameters have been fitted to reproduce hadronic masses and saturation properties 
of groundstate nuclear matter. Within the same model and with the same parameters also 
nuclei and hypernuclei have been described quite successfully. The resulting star shapes 
show a significant excentricity which can also be seen from its hypermatter distribution in 
the interior of the star. The moment of inertia shows a rather strong dependence on the 
rotational frequency, which has to be taken into account for analyzing the spinning-down 
and related age-estimates of pulsars. No back-bending of the inertia occurs due to the lack 
of a phase transition of the chiral model for high densities and low temperatures. In general 
the comparison with a quite different model that contains hyperons shows the robustness of 
the results for the rotating star as well as the relative insensitivity of the star properties on 
very different models with quite distinct hyperon content 0. In this context a discussion of 
the possible effects of populating higher baryonic resonances in the neutron star should be 
performed. Those additional degrees of freedom can change the high-density behavior of the 
system significantly, generating a different structure of the time-evolution of the rotating 
star. Such a calculation might rule in or out specific theoretical descriptions of hadronic 
resonances by studying their effect on neutron star properties. Work in this direction is in 
progress [|18| . 
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